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The unusual X-ray emission of the short Swift GRB 090515: 
Evidence for the formation of a magnetar? 
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ABSTRACT 

The majority of short gamma-ray bursts (SGRBs) are thought to originate from the 
merger of compact binary systems collapsing directly to form a black hole. However, it has 
been proposed that both SGRBs and long gamma-ray bursts (LGRBs) may, on rare occasions, 
form an unstable millisecond pulsar (magnetar) prior to final collapse. GRB 090515, detected 
by the Swift satellite was extremely short, with a Tgo of 0.036±0.016 s, and had a very low 
fluence of 2 x 10~^ erg cm~^ and faint optical afterglow. Despite this, the 0.3 - 10 keV flux 
in the first 200 s was the highest observed for a SGRB by the Swift X-ray Telescope (XRT). 
The X-ray light curve showed an unusual plateau and steep decay, becoming undetectable 
after ^500 s. This behaviour is similar to that observed in some long bursts proposed to have 
magnetars contributing to their emission. 

In this paper, we present the Swift observations of GRB 090515 and compare it to other 
gamma-ray bursts (GRBs) in the Swift sample. Additionally, we present optical observations 
from Gemini, which detected an afterglow of magnitude 26.4 ± 0.1 at T+ 1.7 hours after 
the burst. We discuss potential causes of the unusual 0.3 - 10 keV emission and suggest it 
might be energy injection from an unstable millisecond pulsar. Using the duration and flux of 
the plateau of GRB 090515, we place constraints on the millisecond pulsar spin period and 
magnetic field. 
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1 INTRODUCTION provided by X-ray afterglows, it was possible to identify the opti- 

^ , ,. ^ , cal afterglow and the host galaxies of many GRBs. 

Thirty years after t he discovery of Gamma -ray Bursts (GRBs) by 

the Vela satellites (iKlebesadel e taD ll973h . the first X-ray after- With the launch of the Swift satellite (iGehrels et al.ll2004l) . the 

glow was d etected for GRB 970228 by the Beppo-SAX satellite X-ray afterglow has been studied in great detail, placing fighter 

(ICostaetal.i ri997,) . With the increased accuracy for the position constraints on models for GRB emission. Additionally, Swift has 

enabled the detection of Short GRB (SGRB) X-ray afterglows, al- 
lowing t hem to be dire ctly compare d to Long GRB (LGRB) af- 
terglows JGehrels et al.ll200 5 ; B arth elmv et a . l2005al : iHiorth et al] 



E-mail: bar7@star.le.ac.uk 12005 ah . iNousek et all (120061) and lZhang et all J2006h described the 
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"canonical GRB light curve" as three stages comprising a steep de- 
cline followed by a shallower decay and then a final decay phase. 
lO'Brien et alJ (1200 6) showed that not all X-ray lig ht curves for 
GRBs are of the "cannonical" variety. They and Wil lingale et al.l 
(120070 suggested that the X-ray light curve comprises two main 
components, the prompt emission and the afterglow. The rela- 
tive strength of these components determines the observed X-ray 
light curve. A more recent study of all Swift X-ray afterglows by 
lEvans e t al] (12009) has reinforced these findings. The initial steep 
decay following the prompt emission typically has a power law de- 
cay with index a ^-^ 2 — 5, where / oc t~^ (tis the time after the 
burst in seconds and f is the flux) J O'Brien et al.l l2006). 

Multi- wavelength observations have associated LGRBs with 
type Ibc core coll a pse supernoyae at cosmological distances (e.g. 
lHiorthetal.ll2003l : lstaneketal.ll2003l), although not al l such su- 
pernovae produce long GRBs JWooslev & Bloomll2006h . The pro- 
genitors of SGRBs are less well understood, but the most pop- 
ular theory is that they originate from the merger of compact 
binary systems, for example ne utron stars or a neutron star 
and a black hole J Lattimer & Sch rammlll97^ : lEichleretal.ll 19891 : 
iNaravan. Paczvnski, & Piran 1992). It has also been suggested that 
both LGRB and SGRB progenitors could produce an unstable mil- 
lisecond pulsar. This is expected to contribute a small fraction 
of the GR B population (|Usovlll992l: IPuncan & ThompsonI 19921 : 
lDai&Lul[T998 a.b: Zh ang & MeszarosI 1200 ih . iTroia et al.r(l2007h 
and lLvons et al.l ((2010) found examples of LGRBs that have an ob- 
servable plateau and steep decay in the X-ray light curve, which 
have been intepreted as caused by energy injection by an un- 
stable millisecond pulsar which then collapses. Magnetar models 
have also been proposed to explain late central engine activity in 
SGRBs, for exaniple late tim e plateaus in the X-ray afterglows 
(iFan & Xull2006l : IYu. Chei^. & Cao 2010: iiDair Osso et al. | l201Qh 
and X-ray flares dPan. Zhang. & Progal l2005l : iGao & FanI l2006h . 
Here we present an analysis of GRB 090515 which is the best case 
for an early X-ray plateau in an SGRB. 

GRB 090515 was one of the shortest GRBs observed by Swift, 
with among the lowest fluence, yet for ^200 s it had the brightest 
SGRB X-ray afterglow and did not appear to be fading until a sud- 
den steep decline at ^ 240 s. After the first orbit, it was not detected 
again. Explaining this unusual X-ray behaviour is the subject of 
this paper. We describe the observations of GRB 090515 in section 
2, compare it to other GRBs in section 3 and discuss the poten- 
tial origin of the unusual X-ray emission in section 4. Throughout 
the paper we adopt a cosmology with Hq — 71 kms~^ Mpc~^, 
Q^rn = 0.27, r^A = 0.73. Errors are quoted at 90% confidence for 
X-ray data and at la for optical data. 



2 OBSERVATIONS 

2.1 Swift Observations 

All analysis has been performed by using standard routines in 
HEASOFT, XSPEC, Q DP and the automatic X-ray Telescope (XRT, 
iBurrows et al.ll2005[) data products produc ed by the UK Swift Sci- 
ence Data Centre (lEvans et al.ll2007lJ2009h . 

Swift triggered on GRB 090515 at 04:45:09 UT on 15th May 
2009 with BAT position RA = lOh 56m 41s and Dec = +W 11' 
lH' (iBeardmore eTaP l2009h . The Ultra- Violet and Optical Tele- 
scope (UVOT) enhanced refined XRT position was RA = lOh 56m 
36.11s and Dec = + 14° 26' 30.3'' with an uncertainty of 2.7" 
(lOsborne et al.ll2009l) . 
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Figure 1. The combined light curve for GRB 090515, in grey are the BAT 
data and black are the XRT data. In the lower box there is the hardness ratio 
for the BAT data ((50 - 100) keV/(25 - 50) keV) in grey and the hardness 
ratio for the XRT data ((1.5 - 10) keV/(0.3 - 1.5) keV) in black. Inset is the 
BAT count rate per detector light curve with linear time. 



The Tqo duration of GRB 090515 was 0.036 =b 0.016 s 
(Barth elmv et al.ll2009 a'). The spectrum of the prompt gamma-ray 
emission can be fit by a single po wer law, of photon index F^ = 
1.6±0.2 (JBarthelmv et al.l2009ah . The fluence is 2.0=b0.8 x 10"^ 
erg cm ~^ and the peak photon flux is 5.7 ± 0.9 ph cm~^ s~^. All 
values are in the 15-150 keV energy band. The BAT light curve is 
shown in Figure[T]as the grey data points and also shown in the inset 
with linear time. The BAT count rates were converted to flux in the 
energy band 0.3-10 keV using the average spectral index for the 
BAT and the XRT spectra. There is no evidence of extended emis- 
sion d etected in the BAT energy range (JNorris. Gehrels, & Scargld 
|20i3). 

We completed a spectral lag analysis for GRB 090515 using 
the cr oss correlation function method described in lUkwatta et al.l 
(120100 , the 8 ms time binned lightcurve and BAT channels 1, 2 
and 3. Channel 4 did not detect enough emission to make a lag 
measurement. The lag times are (with la errors): lag(Ch2-Chl)= 
6 ± 4 ms, lag(Ch3-Ch2)= 3 di 2 ms and lag(C h2-Chl)= 10 di 4 
ms. Typically SGRB s have negligble lag times (INorris & Bonnelll 
I2OO6I : lYietal.1 12006^ and LGRBs have typical lag times ranging 
from 20 ms to ^1000 ms (Ukwatta et al. 2010), so it is interesting 
that GRB 090515 appears to have a small lag time. 

The X-ray spectrum in the 0.3 - 10 keV energy band is best 
fit by an absorbed power law with Vx — 1.88 =b 0.14 and A^h = 
6.115° X 10^° cm-^ in excess of the Galactic A^h = 1.9 x 10^° 
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cm~^ (IBeardinore & Evanj|2009h . The X-ray light curve is best fit 
by a broken power law with 2 breaks giving a reduced x^ of 0-86. 
The initial decay is relatively flat (ai = 0.291q§7) with a break at 
Ti = 156.2l^g^2S followed by a steeper decay of a2 = 2.51+070. 
At T2 = 240.8^9 8 s it breaks to an extremely steep decay of 
0^3 > 9. Although, we have fitted the X-ray light curve using a 
broken power law, we note that the decay appears to be a smooth 
curve. The X-ray light curve is shown in Figure[T]as the black data 
points and the lower panel shows the hardness ratio for the gamma- 
ray emission (in grey), i.e. the ratio of the 50 - 100 keV emission 
to the 25 - 50 keV emission, and the hardness ratio of the X-ray 
emission (in black, (1.5-10) keV/(0.3 - 1.5) keV). The hardness 
ratio is fairly constant during the plateau, with the exception of a 
point at ^120 s that could be a flare and does correspond to a small 
peak in the X-ray light curve, but this may just be noise. There 
are insufficient counts to characterise the hardness ratio during the 
decay. 



2.2 Early Optical Observations 

The field of GRB 090515 was observed at early times by several 
optical telescopes but none detected an optical afterglow. The up- 
per limits of the R band and white filter observations are given in 
Table 1 . During the plateau phase, we can predict the optical flux 
density, assuming that the X-ray and optical emission are from the 
same emitting region. If there is not a cooling break in the spec- 
trum (i.e. Tx — Tox) then we would expect the optical flux to 
be 1.7 X 10~^^ erg cm~^ s~^ Hz~^, corresponding to an appar- 
ent magnitude of i? = 15.6. This is brighter than all of the op- 
tical upper limits during the plateau, so we should have observed 
the optical afterglow. However, if there were a cooling break in 
the spectrum between optical and X-ray then Vox — ^x — 0.5 
and, in this case, the optical flux density would be 8.7 x 10~^^ erg 
cm~^ s~^ keV~^Hz~^, corresponding to an apparent magnitude 
of R — 21.3. This is slightly deeper than the optical upper limit 
provided by UVOT. Therefore, if the optical emission was from the 
same emitting region as the X-ray and there is a cooling break in 
the spectrum, there is a slim chance that the optical flux was be- 
low the observed limits so the non-detection is consistent with the 
X-ray data. 



2.3 Gemini Observations 

We obtained optical observations of GRB 090515 using Gemini 
North and GMOS beginning at 06:26 UT, approximately 1.7 hours 
after the burst, with a second epoch observation being taken on the 
subsequent night, and a final comparison epoch on 28 November 
2009. The images were obtained in the r-band , and were reduced 
via the standard IRAF Gemini tasks (lTodvlll993h . The image condi- 
tions for our first epoch were excellent, with seeing of 0.5'^ result- 
ing in extremely deep imaging in our total exposure time of 1800s. 
A full log of observations is shown in Table [21 

Within the refined XRT error circle we locate a single source at 
RA = lOh 56m 35.89s and Dec = -hl4° 26' 30.0'', with a magnitude 
of r = 26.36 =b 0.12, calibrated against existing SDSS observa- 
tions of the field, shown in Figure [2l This source is still visible, but 
at lower significance in our shallower images obtained on 16 May 
(r = 26.54 zb 0.33). In our final epoch there is no source visible 
at the afterglow location, to a limiting magnitude of r > 27.4 con- 
firming a fading counterpart. We therefore conclude that this is the 
optical afterglow of GRB 090515. At r = 26.36, this is the faintest 
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Figure 2. The circle marks the location of the XRT enhanced position of 
GRB 090515 on the Gemini observation from epoch 1 on the left and epoch 
3 on the right. An optical afterglow candidate is observed within the error 
circle. Labeled are the brightest nearby galaxies. 



Table 3. Photometry of the nearby galaxies, as labeled in Figure [21 cali- 
brated using SDSS observations. 



Galaxy 


Magnitude 


Offset (") 


1 


20.2 ±0.1 


14 


2 


21.3 ±0.1 


16 


3 


22.5 ±0.1 


15 


4 


22.6 ±0.1 


13 


5 


23.4 ±0.1 


6 



GRB afterglow ever discovered at similar times after the burst, and 
confirms the necessity of rapid and deep observations with 8-metre 
class observatories. As the observed X-ray absorption is relatively 
low (A^H '^ 6 X 10^° cm~^), the faint optical afterglow is unlikely 
to be a consequence of extinction (unless it is at high redshift). The 
optical afterglow has a relatively flat lightcurve, with a decay slope 

ofo.o6l°;?9- 

Comparing this afte rg low to the sample in 
iNvsewander. Fruchter. & Pe'ed (120091) . we note that this is 
the first SGRB with a fluence below 10~^ erg cm~^ with a 
detected optical afterglow. Additionally, the afterglow at 1.7 hours 
is fainter than all the observed optical afterglows at 1 1 hours. GRB 
080503 also had an initially very faint optical afterglow, but it then 
rebrightened t o a peak of r ^ 25 .5 at 1 day and no host galaxy 
was identified (' Perlev et al.ll2009ah . 

Assuming there is not a cooling break in the spectrum, i.e. 
Tx = Tox, we predict that the X-ray flux, 0.3 - 10 keV, at the 
time of the optical observations should be 6.6 x 10~^^ erg cm~^ 
s~^. This is consistent with the observed upper limit. 

Labeled in Figure |2] are the five brightest nearby galaxies and 
Table [3] provides their magnitudes and offsets from the GRB lo- 
cation. These galaxies are candidates for the host galaxy of GRB 
090515, with significant offsets, or the burst could be associated 
with a significantly fainter underlying host galaxy. 



3 COMPARISON TO OTHER GRBS 

The XRT fight curve of the low fluence GRB 090515 is unusual 
as it goes from being the brightest SGRB in X-rays to one of the 
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Table 1. The optical observations of the field of GRB 090515. 



Telescope Mid point time after trigger 

(s) 



Exposure Time Band 



Upper Limit 
(magnitude) 



Flux Upper Limit 

(erg cm~^s~-'^Hz~-'^) 



KAIT 


20 


Super LOTIS 


43 


ROTSE III 


86 


UVOT 


142 


UVOT 


1228 


KAIT 


2078 


Lick 


2286 


ROVOR 


5496 



540 


R 


19.1(1) 


10 


R 


17.7(2) 


67 


R 


18.4(3) 


146 


White 


20.35(4) 


488 


White 


21.24(4) 


540 


R 


20.5(1) 


60 


R 


21.3(5) 


4200 


R 


21.4(6) 



6.6 X 10-28 
2.4 X 10-2'^ 
1.3 X 10-2'^ 

2.1 X 10-28 

9.2 X 10-29 

1.8 X 10-28 

8.7 X 10-29 

7.9 X 10-29 



(i) |Lietal.l ( Eooi) . (2) Iwilliams et al.' ('2009'), (3) 'RuioDakarn et al' ('2009'), (4) 'Seigel & Beardmord ( Eo09h . 

(5)Iperley, Kisla k & Qaneshaling am i 2Q09b ,), (6) Paceetal.. (2009^). 



Table 2. Log of Gemini observations. 



Epoch 


Date start 
(UT) 


Start time after trigger 

(s) 


Exposure time 

(s) 


Filter 


Seeing 
(arcsec) 


Airmass 


Magnitude 


Flux 

(ergcm-2s-iHz-i) 


1 

2 
3 


May 15 06:27 UT 

May 16 05:44 UT 

November 28 14:20 UT 


-6100 
-9x10^ 
- 1.6x10'^ 


1800 
1800 
2800 


r 
r 
r 


0.5 

1 

0.8 


1.021 
1.005 
1.226 


26.36 ±0.12 

26.54 ± 0.33 

>27.4 


8.2 X 10-31 

6.95 X 10-31 

< 4.55 X 10-31 



Table 4. The GRBs considered in detail in this paper. 



10-9 ^ 



GRB 


T90 
(s) 


r(15-150keV) 


Fluence(15-150keV) 
(10-8 erg cm-2) 


090515 


0.036±0.016 


1.6±0.2 


2.0±0.8(i)'(2) 



090607 

080520A 
080503 

070724A 
070616 
070209 
060717 

051221B 
051105 
050813 

050509B 
050421 



2.3±0.1 

2.8±0.7 

170±20 

0.4±0.04 

402±10 

0.10±0.02 

3.0±1 

61±1 

0.028±0.004 

0.6±0.1 

0.048±0.022 

10.3±2 



1.25±0.30 
2.90±0.51 
2.00±0.13 
1.81±0.33 
1.61ib0.04 
1.55±0.39 
1.72±0.38 
1.48±0.18 
1.33±0.35 
1.19±0.33 
1.5±0.4 
1.7±0.4 



11±2(3) 
5.5±1.4(4) 
200±10(5) 
3.0±0.7 (6) 
1920±30 C^) 
1.1±0.3 (8) 
6.5±1.6(9) 
113±13(io) 
2.0ib0.46 (11) 
4.4=bl.l (12) 
0.78±0.22 (13) 
8.8±2.9 (14) 



(1) iBarthelmv etli] (|2009ah (2) I S akamoto & Beardmord (l2009h (3) 
'Barthelmv et al. {2 mm~^ ISakamoto et al.l (l2008h (5) lUkwatta et al.l 
12008.) (6) .Parsons et all (t2007l) C^) | Sato & Barthelmv! (l2007h (8) 

enimore 
120051) 



ISakamoto et al.N2007l) (9 ) Markwar lTet al.l ( I2OO6I) (1°) iFenimore et al.1 
(l2005h (11) iBarbier et al. (.20051) (12) ISatoetal.1 



iBarthelmv et al.l ( [2005bl) (i4) |Sakamoto et"ZI 120051) 



faintest within seconds. The fluence in X-rays during the plateau is 
significantly higher than the fluence in gamma-rays. Additionally, 
the final decay is the steepest decay observed to date (IE vans et alJ 
12009.) . The X-ray spectral index of GRB 090515 is not unusual 
compared to other SGRBs. In Tabled we provide a summary of 
the properties of the long and short GRBs to which we compare 
GRB 090515 in detail 

In Figure [21 we show the 15 - 150 keV fluence and 0.3 - 10 
keV flux at to + 100 s for all the SGRBs in the Swift sample with 
T90 ^ 2s and which were observed by XRT at this time. GRB 
090515 is shown with a filled circle. As expected, the higher fluence 




10^ 



10^ 
15 - 1 50 keV Fluence (erg cm-^) 



1(^^ 



Figure 3. The fluence in the energy band 0.3 - 10 keV versus the 15 - 150 
keV flux for all Swift SGRBs which were observed at 100s after the trigger 
time. The filled circle marks the location of GRB 090515. 



GRBs tend to have higher flux X-ray afterglows. GRB 090515 is 
an exception to this alongside GRB 070724A; both of these bursts 
have an unusually high initial X-ray flux for their fluence. In Fig- 
ure[5la), we compare the combined BAT-XRT light curves of GRB 
090515 and GRB 070724 A. The initial XRT flux of 070724A ap- 
pears to be consistent with flares (as there is a varying hardness 
ratio) and an underlying broken power law decay. There is no ob- 
vious plateau phase for GRB 070724A, but this may have occurred 
prior to the XRT observations. The steep decay phase of GRB 
070724A, with a = 3.44+0-35 is much shallower than the steep 
decay of GRB 090515. Additionally, the optical afterglow of GRB 
070724A had a magnitude of i = 23.79=b0.07 at 2.3 hours after the 
burst, corresponding to a flux o f 6.86 x 10~^° ergcm~^ s~^ Hz~^, 
with an associated host galaxy (JBerger et al.ll2009al : lKocevski et al.l 
120101) . This flux is almost an order of magnitude larger than the op- 
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Figure 4. (a) The optical flux light curves for all observed SGRB af- 
terglows in the R band, (b) Normalised using the XRT flux at 1000 s. 
Co lour scheme : GRB 090515 - purple star, GRB 091109B - smaU green 
X (lLevanetal.1 2009: Malesani et al. 2009), GRB 090426 - dark blue open 
square jAntonelli et al. 2009; Xin et al. 2010), GRB 090305 - light blue 
open diamond Jcenko et al. 2009; Berger & Kelson 2009b), GRB 080905A 
- small red open circle ( Rowlinson et al. 2010), GRB 071227 - green 
fill ed triangle j Berger, MorreU, & Roth 2007), GRB 070809 - large black 
X dPerlev. Thoe ne. & Bloom 2007; Perlev et al. 2008), GRB 061201 - large 

S )ink open cir cle ( Stratta et al. 2007), GRB 060121 - dark grey fill ed circle 
Levan et al. 2006), GRB 05 1221 A - orange open cross ( Soderb erg et all 
l2Q06h . GRB 050709 - fight grey filled square iHiorth et al.„20Q5a) . (c) The 
optical flux fight curve for GRB 090515 (purple stars) with GRB 080503 
(red). 



tical afterglow of GRB 090515 at 1.7 hours and GRB 090515 does 
not have an identified host galaxy. However, GRB 070724A does 
share many similarities with GRB 090515 so we cannot rule out 
the possibility that they originate from a similar progenitor. 

Figure l^a) shows the lightcurves for the observed R band 
optical afterglows associated with SGRBs (published values con- 
verted from magnitudes into flux density in Jy), GRB 090515 is the 
faintest observed and one of the earliest detections after the trigger 
time. In Figure |4tb) we have divided the optical fluxes by the XRT 
flux at 1000 s after the trigger time. When we have considered the 
XRT flux at 1000 s, the optical afterglow of GRB 090515 is not un- 
usually faint compared to other SGRBs. We also show the optical 
light cu rve for GRB 080503 (a short burst with extended emission 
IPerlev et al...2009a) in Figuregfc) in comparison to GRB 090515. 



3.1 GRBs with similar fluence to GRB 090515 

As the fluence of GRB 0905 1 5 in the 1 5 - 1 50 keV energy band was 
one of the lowest fluences observed for SGRBs, here we compare 
it to other low fluence GRBs. 

GRB 050509B and GRB 050813 were short GRBs detected 
by the Swift satellite that were similar to GRB 090515 during the 
prompt emission phase. However, the combined BAT and XRT 
light curves for GRBs 050509B and 050813, shown in Figure 
[5lb), do not show the same X-ray plateau extending to ^200 
s after the burst. GRBs 050509B and 050813 have both been 
used to place constraints on the compact binary merg er model 
of SGRBs (Gehrels et al. 2005; Hjorthetal. 2005b; Bl oom et all 
l2006l : iFerrero et al.1 l2007h . The observed upper limits for GRB 
090515 at late times (after 400s) are consistent with the later emis- 
sion observed for GRBs 050509B and 050813. This suggests that 
the plateau and steep decay are an additional component in the light 
curve of GRB 090515. 

GRB 051105 is a SGRB with an identical fluence to GRB 
090515, but its afterglow w as undetectable by X RT in observations 
starting 68 s after the burst (iMineo et al.ll2005a ll. GRB 070209 had 
the lowest SGRB fluence and was also undetectable by XRT in 
observations starting 78 s after the burst (Sato et al. 2007) . 

In Figure [SJc), the X-ray light curve of GRB 090515 is 
compared to the two lowest fluence LGRBs in the Swift sample 
which were detected by XRT. These are GRB 080520A and GRB 
0607 17 A, they both have significantly higher fluence in the 15 - 
150 keV band than GRB 090515 (due to having longer durations), 
but are a lot fainter in X-rays, again suggesting additional X-ray 
emission in GRB 090515. 

It is possible that these GRBs had plateau phases which 
end prior to the XRT observations. However, as Swift slewed 
promptly to these GRBs (observations typically starting within 100 
s), a plateau phase would need to be significantly shorter than 
that observed for GRB 090515. The main exception to this is 
GRB 0607 17 A, which had XRT observations begining when GRB 
090515 was in the steep decay phase. 



3.2 GRBs with steep decays 

The GRB with an X-ray light curve most similar to GRB 090515 is 
GRB 090607, which has a Tgo just above the short-long boundary. 
They are compared in Figure|5ld). Both light curves show a distinc- 
tive steep decay at ^200s. However, the emission of GRB 090607 
between 80 and 100 s is not a plateau as observed in GRB 090515 
and, given the hard spectrum which softens as the emission decays 
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Figure 5. The BAT-XRT light curve and hardness ratios for GRB 090515 in blue in comparison to other GRBs. (a) GRB 070724 A in red. (b) ORB 050813 in 
red and GRB 050509B in purple, (c) GRB 060717A in red and GRB 080520A in purple, (d) GRB 090607 in red. (e) GRB 050421 in red and GRB 080503 in 
purple, (f) GRB 070616 in red. In the lower boxes for each graph, there is the hardness ratio for the BAT data ((50 - 100) keV/(25 - 50) keV), with a star, and 
the hardness ratio for the XRT data ((1.5 - 10) keV/(0.3 - 1.5) keV). 
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(shown in the lower panel of Figure [3d)), is more likely due to a 
flare at the start of the XRT observations. 



4 DISCUSSION 

The steep decay in the unusual X-ray light curve of GRB 090515 
cannot be explained using the external shock afterglow models. In- 
stead, we consider if this GRB was a naked burst with faint, rapidly 
fading emission, and if the X-ray plateau is powered by an unstable 
millisecond pulsar. These possibilities are discussed below. 



4.1 An under-luminous naked LGRB 

If a GRB occurs in a very low density ISM then the afterglow 
from external shocks between the jet and the ISM could be too 
faint for detection by Swift. Instead, there would just be the prompt 
emission followed by a rap id decline due to the "curvature effect" 
(iKumar & Panaitescul200Qh . This predicts a decay in flux described 
by: 



/. oc v-^t-^-^ 



(1) 



where /3 is the observed spectral energy index at frequency z^ (/3 = 
r — 1), and t is the time since the trigger. We should observe a decay 
a = 2 + /3 = r+l.We compare here GR B 090515 with a good 
candidate for a naked burst, GRB 050421 (Godet et al.ll2006h . 

GRB 050421 was a weak long GRB detected by BAT fol- 
lowing a steep decay, as shown in Figure Oe), although the de- 
cay is not as steep as for GRB 090515. There is evidence of spec- 
tral evolution, as the emission is getting softer (the lower panel of 
Figure [3e)); however, the spectral evolution is during the steep 
decay and not the plateau region. The initial hardness ratio for 
(1.5 - 10) keV/(0.3 - 1.5) keV i s 6 tim es larger for GRB 050421 
than GRB 090515. iGodet et all ([200^ explained the steep decay 
(a = 3.1 ± 0.1) of GRB 050421 by assuming it was a "naked 
burst", i.e. there was no forward shock component of the after- 
glow as the interstellar medium (ISM) was not dense enough for the 
shock wave to produce a typical afterglow. The detected decaying 
emission is consistent with the "curvature effect". GRB 080503, 
shown in Figure [5je), has also been explained as a short "naked 
burst" with extended emission detected in the BAT (although not a 

S >lateau), the X-ray decay is consistent with the "curvature effect" 
Perlev et al.l2009al) . However, the steep decay for both of these are 
significantly shallower than the decay of GRB 090515, which was 
Q^3 >9 (with to at the start of the prompt emission). 

GRB 090 515 shares some similarities with GRB 050421 and 
GRB 080503 (iGodet et al.ll2"00^ : IPerlev et al.ll2009ah . lZhang et all 
(l2009h suggested that the burst duration, observed by BAT, rep- 
resents the duration that the jet is relativistic and, with a non- 
relativistic (or less relativistic) jet, the central engine can be active 
for longer than this time and may be observed by XRT. Therefore, 
the X-ray plateau observed for GRB 090515 could be a continua- 
tion of the prompt emission, which has fallen below the threshold 
of BAT. So with a more sensitive detector, GRB 090515 may have 
been identified as a LGRB. If true, we should expect to see that 
the steep decay matches the "curvature effect" like GRB 050421. 
During the plateau, the spectral index Tx is 1.88 ± 0.14 predicting 
a steep decay slope of a — 2.88 =b 0.14. As the observed decay 
is significantly steeper than this, it does n ot fit the "curvatur e ef- 
fect" theory. Using the method described bv lLiang et al .*bOO^, we 
shifted the to to the possible flare at the end of the plateau in GRB 
090515. The steep decay becomes less extreme, a = 3.7 ± 0.6, 



but still steeper than the predicted decay slope. This method relies 
on correctly identifying the time at which the central engine is last 
active and with a plateau in the light curve this point is difficult 
to identify. The steep decay of GRB 090515 following the plateau 
may be consistent with the "curvature effect" if a later location of to 
is identified. Alternatively, this could be associated with a narrow 
opening angle for the jet which creates the plateau, as in that case 
outside of l/F there would be very little high latitude emission, giv- 
ing a much steeper decay slope. It is also possible that the spectrum 
softens immediately prior to the steep decay, however we do not 
have enough observed counts to produce a reliable X-ray spectrum 
at this time. 

GRB 090515 can potentially be explained as an under- 
luminous naked long GRB, however this is reliant on the assump- 
tion that the plateau is powered by prolonged activity in the central 
engine. 



4.2 An unstable millisecond pulsar (magnetar) central engine 

The bright X-ray plateau in the light curve of GRB 090515 
could be associated with the formation, emission and collapse 
of a millisecond pulsar. There have been predictions that in 
some GRBs an unstable millisecond pu ls ar may be forni ed 
(Usov 1992; Dunca n_& Thom pson 1992; Dai & Lul Il998al lbl: 
Zhang & Meszaros 200l|). At formation, there is enough rotational 
energy to prevent gravitational collapse. This energy can be re- 
leased as electromagnetic radiation or gravitational waves, causing 
the pulsar to spin down until it reaches a critical point at which it is 
no longer able to support itself. At this point the pulsar collapses to 
a black hole and the emission stops. This would be evident in the X- 
ray light curve as a plateau caused by energy injection from the mil- 
lisecond pulsar followed by an extremely steep decay when the pul- 
sar collapses. We might expect millisecond pulsars formed during 
the core collapse of a massive progenitor st ar to be associated with 
long GRBs and th is has been suggested bv lTroiaet al.l (120071) and 
lLvonsetal.l ( l201Ql) . GRB 090515 was an extremely short GRB, but 
a millisecond pulsar could be formed by two merging neutron stars 
(a potential progenitor of SGRBs), depending o n various assump- 
tions about the neutron stars' equa tions of state JPai & Lulll998al : 
iDai et al.ll2006l:[Yu & Huandl2007h. 

iTroia et all ( 120071) and lLvons et al.l (|201Q|) studied LGRBs 
with a plateau and a steep decay and GRB 090515 shows simi- 
larities to them. In Figure Qf), we compare the light c urve of GRB 
090515 to that of GRB 070616 (IStarling et all 120081) . one of the 
sample chosen by iLvons et al.l (120101) as potentially showing evi- 
dence of an unstable millisecond pulsar. When comparing the light 
curves, GRB 070616 appears to be a brighter and longer version of 
GRB 090515 but with a bright afterglow component at later times. 

We have used the following equations from 
IZhang & MeszarosI (120011) (equations [2 and [3]) to determine 
if GRB 090515 could be a millisecond pulsar, using Tem,3, the 
rest frame duration of the plateau in units of 10^ s, and Lem,49, 
the luminosity of the plateau in units of 10^^ erg s~^, in the rest 
frame energy band 1 - lOOOkeV. The equations are rearranged to 
give equations |4] and [5] these are used to predict the magnetic field 
strength and the spin period of a pulsar formed by this method. 



Te^,3 = 2.05 (/45B-?5Po'-3i?6"') 

i^em,49 ^ {Bp^iqPq_^Rq) 

^p,15 — 4.2025/45/^6 /^em,49^em,3 

^0,-3 — 2.05/45/^em,49^em,3 



(2) 
(3) 
(4) 
(5) 
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where /45 is the moment of inertia in units of 10^^ g cm^, Bp,i5 
is the magnetic field strength at the poles in units of 10^^ C, Rq 
is the radius of the neutron star in 10^ cm and Po,-3 is the ini- 
tial period of the compact object in milliseconds. These equations 
apply to the electromagnetic dominated spin down regime, as the 
gravitational wave dominated regime would be extremely rapid and 
produce a negligble effect in our analysis. We could assume, as in 
iLvons et aLI (120101) . that we can use standard values for a neutron 
star so that Ia^ ^ 1 and Rq ^ 1 which may be appropriate for a 
collapsar. However, as we would be forming an unstable millisec- 
ond pulsar by merging two neutron stars the true values may be 
different, depending on the mass and equation of state. For a mil- 
lisecond pulsar formed by a binary merger, we can take the mass 
of the neutron star to be Mns =2.1M0 (iNice et all l2005 l) and 
estimate /45 '^ 1-5. Although GRB 090515 has many properties 
similar to other SGRBs suggesting the progenitor is most likely a 
compact binary merger, there have been predictions that collapsars 
may also produce a SORB (for examp le from an orphan precursor 
jet, Uaniuk. Moderski. & Progall2008h and evidence that a signifi- 
cant fraction of SGR Bs are related to collapsars rather than com- 
pact binary mergers ("Virgili et al]|2009l : ICui. Aoi. Nag ataki, & Xul 
I2OI O). So in the following analysis we compare both progenitor 
models. 

As a redshift was not obtained for this GRB, we used a range 
of redshifts from z — 0.2 up to an upper limit of z — 5.0 con- 
sistent with the detection of the optical afterglow. We assume that 
the millisecond pulsar was formed at t ^ and, hence, the dura- 
tion of the plateau in the observer frame is 240 s. We calculate the 
luminosity of the plateau using the observed 0.3 - 10 keV flux of 
^ 1 X 10~^ erg cm~^ s~ \ the spectral index during the plateau 
(1.88) and a k-correction (B loom. Frai L & Sari 2001). These val- 
ues were then substituted into the equations ^ and (|5]) to calculate 
Bp,i5 and Po,-3. These are plotted as a blue contour in Figure[6ta) 
assuming it was formed from a collapsar and a purple contour if 
formed by a binary neutron star merger. 

Also shown in Figure[6{a), are the region s in which a millisec - 
ond pulsar would be expected, as defined in iLvons et al.l (I2OIC1I) : 
the red line represents the breakup spin-period for a neutron star of 
mass I AM(^ (^0.9 6ms, 'Lattimer & Prakashll2004h . Using equa- 
tion [6l ( Lattimer & Prakash 2004), we calculate this limit for the 
binary merger scenario with a mass of 2.1 Mq to be P^O.66 ms 
(where P is the minimum spin period of the neutron star in ms) and 
this is shown with a red dashed line. 



^ 0.81M^^^"R\ 



l/2p3/2^ 



(6) 



The initial rotation period needs to be ^lOms ( Usov 19921) an d 
the lower limit for the magnetic field is ^ lO^^G (! Thompsonll2007n . 
This shows that GRB 090515 could have formed a millisecond pul- 
sar if it had a redshift of 0.3 < z < 3.5 for a collapsar progenitor or 
a redshift of 0.2 < z < 4.4 for a binary merger progenitor. These 
are both very reasonable redshift ranges when we compare them 
to the sample of GRBs. The magnetic field for a given spin period 
is slightly lower for a binary merger progenitor than for a collap- 
sar progenitor. Alo ngside the prediction by Troia et al.l (l2007h and 
iLvons et al.l (120 101) of a plateau followed by a steep decay for the 
lightcurve of a millisecond pulsar collapsing to a black hole, which 
matches the observed light curve for GRB 090515, this analysis 
provides a consistent case for GRB 090515 forming a millisecond 
pulsar irrespective of the two initial progenitor models considered. 
Using a causality argument, i.e. that the speed of sound on the 
neutron star cannot exceed the speed of light, we can place a tighter 
constraint on the minimum possible radius, Rq and Mi. 4, where 




Period (ms) 



Figure 6. (a) The blue line shows the magnetic field and period for a mil- 
lisecond pulsar formed during GRB 090515 as a function of redshift assum- 
ing a neutron star mass of 1.4M0 and the purple dotted line assumes a neu- 
tron star massof24Mo. The g reen stars are the 18 degree beamed LORE 
sample from lLvons et al.l ( I2OIOI) . The red line shows the limit at which the 
progenitor would violate the breakup spin period of a pulsar for a mass of 
1.4M0 and the dashed red line is f or a mass of 2.1M0. The other regions 
are as defined in 'Lyons et al.l ( |201Q|) : dark grey shading corresponds to for- 
bidden regions (assuming a mass of I.4M0) and light grey are limits based 
on the previous studies (as discussed in the text). The dotted lines represent 
contours of equal redshift decreasing from left to right, (b) The upper mag- 
netic field limit in (a) has been replaced by the red curved line giving the 
forbidden region assuming causality shaded in dark grey (assuming a mass 
of I.4M0). This region will change depending on the mass of the neutron 
star, the rest frame duration and luminosity of the plateau. The red dashed 
curved line represents the forbidden region for a binary merger progenitor. 
(c) The different contours represent the effect of increasing the radius of the 
neutron star from 10 km to 30 km assuming a constant mass of 1.4M0. 
Additionally, we include a limit imposed on redshift due to detection of the 
afterglow in the R-band. 
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Ml A is the mass of the neutron star in IAMq, using equation [7] 
(ILattimer et alJI 19901) . The moment of inertia, given in equation [S] 
is based on the assumption that the neutron star can be modelled as 
an uniform sphere. 



Re > 0.6225Mi.4 
/45 - MiaRI 



(J) 
(8) 



This constraint on radius and moment of inertia for a given 
mass can be substituted into equations |4] and [5] to define a forbid- 
den region for a given neutron star mass, plateau duration and lu- 
minosity. The forbidden region is described by equations [9l and [Tol 
and is shown in Figure [S^b) for GRB 090515 assuming a mass of 
1.4M0, for a collapsar progenitor (red curved line), and 2.IM0, 
for a binary merger progenitor (red curved dashed line). 



^p,i5 > lO-STg^ 3 ^ 



'em, 49 
-1 



(9) 
(10) 



It has been suggested that the radii of proto n eutron stars may 
be as large as a few tens of kilometers dOtt et al.|[!200 6), so in Fig- 
ure [6[c) we show the effect of increasing the radius, from 10 km 
to 30 km, for a mass of IAMq, using the plateau luminosities and 
durations previously calculated for GRB 090515 assuming it is at 
a range of redshifts. For larger radii, the unstable millisecond pul- 
sar has to be at higher redshifts, have a smaller magnetic field and 
larger period. As we have an R-band detection of the optical after- 
glow, we can place the upper-limit z^5 on the redshift . 

In Figure [T] we inv estigate the effect o f the different beam- 
ing angles considered by iLvons et al] (I2OIO ) assuming a mass of 
IAMq. As the causality forbidden region shown in Figure[6tb and 
c) also depends on beaming angle w e have reverted to using the re- 
gions defined bv'L vons et al.l (|201Q|) for clarity. Up to this point, we 
have only considered isotropic emission and this shows beaming 
the emission would greatly affect the results obtained. Simulations 
have shown that a relat ivistic jet can be produced by a magnetar 
(IBucciantini et al.ll2009l) . If the emission was beamed by 4 degrees 
the observations would not support the magnetar model. With a 
beaming angle of 18 degrees, GRB 090515 would need a redshift 
of 1 < 2; < 5 in order to satisfy the model and the constraints 
obtained by observing an optical afterglow. The more tightly the 
emission is beamed, the higher the redshift that the burst would 
need to be at in order to fit the magnetar model and this may ex- 
plain why a host galaxy has not been identified. 



5 CONCLUSIONS 

GRB 090515 is a very unusual SGRB, as its low gamma-ray flu- 
ence would lead us to expect a significantly fainter X-ray light 
curve than observed at early times. Most importantly, the X-ray 
plateau followed by an extremely steep decay is very unusual, but 
may not be unique in the Swift sample. With a more sensitive detec- 
tor, the plateau observed by XRT may have instead been identified 
as part of the prompt emission and GRB 0905 15 might instead have 
been classified as a LGRB. Therefore, it poses interesting questions 
about the progenitor model and for the classification of other GRBs. 
In this paper, we have considered the two popular progenitor mod- 
els for GRBs, collapsars and compact binary mergers. 

GRB 090515 is the first SGRB with a fluence below 10"'^ 
erg cm~^ with an observed optical afterglow at 1.75 hours 




10 
Period (ms) 

Figure 7. We show here the effect of beaming the emission of GRB 090515 
assuming a mass of 1.4M0. Solid lines show isotropic solution and th e 
solutions for the two beaming angles considered in Lyons et ^ |2010|) . 
The dotted lines represent contours of equal redshift decreasing from left 
to right. The forbidden regions are as defined for Figure |6la). 



(R=26.4zb0.1), and this is the faintest detected optical afterglow 
for a GRB at that time. 

We suggest that the simplest explanation for the unusual light 
curve of GRB 090515 is that it shows prolonged emission from an 
unstable millisecond pulsar, followed by an extremely steep decay 
when the millisecond pulsar collapses. Given the short duration of 
the GRB and the other properties, we favour the binary merger pro- 
genitor but cannot rule out a collapsar progenitor. For a collapsar 
progenitor, the proposed unstable millisecond pulsar with a spin 
period of 10 ms would have a magnetic field of^ 3 x 10^^ Gatz 
^0.3 and with a spin period of 1 ms the magnetic field would be 
^ 6xl0^^Gatz^3.5. The binary merger progenitor model gives 
a spin period of 10 ms and a magnetic field of^ 2.5 x 10^^ Gatz 
^ 0.2 to a spin period of 66 ms and a magnetic field of'-^4xl0^^ 
G at z ^ 4.4. These values assume isotropic emission and a radius 
of 10 km. 
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